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. ABSTRACT 

0\ We propose a new technique, which we caU the lens paraUax method, to deter- 

mine simuhaneously the redshift distribution of the faint blue galaxies and the 
mass distribution of foreground clusters of galaxies. The method is based on 
O gravitational lensing and uses the fact that the mean angular size of the faint 

2 blue galaxies is a well-determined function of surface brightness. The method 

7—i requires first a calibration of the angular-size vs. surface brightness relation in 

^~~^ unlensed blank fields and a determination of the mean redshift of the brightest of 

the faint blue galaxies. By combining this information with observations of the 
.^ distorted images of background galaxies in the fields of foreground clusters, we 

O show that it is possible to obtain the mean redshift of the galaxies as a function 

I of their surface brightness. With a sample of about ten moderate redshift rich 

■^ clusters and using ten bins in surface brightness, a bin-to-bin signal-to-noise ratio 

^ of ~ 3.5 can be achieved. The method simultaneously allows a determination of 

'^ the convergence and shear of the lensing cluster as a function of position, and 

I through this the mass distribution of the lens can be obtained. The lens parallax 

method can be used in conjunction with, and will improve the accuracy of, other 
"oo previously proposed techniques for mapping clusters. 

ci 

Subject headings: Cosmology: Gravitational Lensing - Galaxies: Clustering - Galaxies: 
Distances and Redshifts - Methods: Observational 

INTRODUCTION 

The effect of gravitational lensing on the images of background galaxies is twofold: the 
images are magnified, and they are sheared. However, quantities like surface brightness 
and colors are preserved. The techniques described in this paper are based on these facts. 
We use a measurement of surface brightness to identify the sub-class to which a particular 
lensed galaxy belongs. We then compare the observed angular size and elliptic distortion 
of the galaxy image to what we expect in the absence of lensing and thereby determine 
the redshift of the galaxy as well as the properties of the lens. 

Observations show that there is a well-defined relation between the angular scale 
length 6s and the apparent B magnitude of the faint blue galaxies, that is to say, the 
scatter around the average angular size 6s{B) at a given B is fairly small, A6s{B) <^ Os{B) 
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(Tyson 1994). Equivalently, there is a tight relation between 6s and the (blue) surface 
brightness /. Gravitational lensing conserves / but multiplies the scale size 6s by ^^/I, 
where /i is the magnification induced by the lens. In addition to an overall magnification, 
the shear in the lens also induces a distortion in the shape of the image. 

For source redshifts Zs ^ 0.5, the surface brightness / of the faint blue galaxies de- 
creases with increasing Zs (e.g., Koo & Kron 1992, Tyson 1994). At the same time, for 
a given lens, the lens-induced magnification and shear of the galaxy images increase with 
increasing Zs- Thus, for a given foreground lens such as a cluster of galaxies, the effects 
of lensing increase monotonically with decreasing /. By measuring the sizes and shapes of 
the faint blue galaxies as functions of / we can thus establish a "redshift ladder" for the 
galaxies. The ladder can be calibrated by direct spectroscopic measurement of the mean 
redshift of the brightest faint blue galaxies. Using this, the mean redshift, Zs{I)^ of galaxies 
at fainter / can then be obtained. This in essence is the lens parallax method which we 
describe in this paper. A virtue of the method is that measurements from many clusters 
can be combined to reduce the noise in the derived Zs{I). 

As a by-product, the method also allows a superior reconstruction of the mass dis- 
tribution of the lensing clusters. Current methods (Kaiser & Squires 1993, Broadhurst, 
Taylor & Peacock 1994, Seitz & Schneider 1994) suffer because they lack information on 
the redshifts of the lensed sources. This problem will be eliminated once the relation 
Zs{I) is obtained using our method. Furthermore, many current techniques (e.g. Kaiser 
& Squires 1993) make use only of the shear-induced distortion of the image and need to 
carry out an integral transform of the observations to calculate the surface mass density 
of the lens. In our method, we measure both the shear and magnification of the lens and 
therefore we obtain direct local measurements of the lens surface density. In this context, 
our method is similar to that described by Broadhurst et al. (1994). 

MEASUREMENTS 

The size and shape of a galaxy image can be quantified by the scale length 6s of its surface- 
brightness profile, the ellipticity x, and the position angle (p^ of the elliptic distortion. We 
adopt the common definition, 

a^ -|- 6^ 

where a and b are the semi-major and semi-minor scale lengths, and we write xi = 
Xcos(2(fi^)^ X2 = X siii(2')^x)- -^^^ ^^^ surface brightness /, we use either the central 
surface brightness or, more conveniently, the total fiux divided by ab. Both quantities are 
invariant under lensing. 

The method proposed here requires that 6s and / be measured first for a large num- 
ber of faint blue galaxies in blank fields. These measurements will yield a relation 6s,i{I) 
between the mean 6s of the field galaxies and their surface brightness /, and also the 
dispersion in values of 6s at a given /. This provides the calibration for later measure- 
ments of galaxies imaged by clusters. It is also desirable that the ellipticity parameter 
X be measured for the field galaxies so as to obtain an estimate of the range of intrinsic 
ellipticities. 



In addition, at the bright end of /, the redshifts of an unbiased sample of faint blue 
galaxies should be measured in order to obtain the mean redshift Zs{I) at this value of /. 
Substantial progress has been made in this direction (Peterson et al. 1986, Broadhurst, Ellis 
& Shanks 1988, Metcalfe et al. 1989, CoUess et al. 1990, Cowie et al. 1990, Guhathakurta, 
Tyson & Majewski 1990, Lilly, Cowie & Gardner 1991, Dressier & Gunn 1992, Koo & 
Kron 1992), and it is possible that the information available now is already sufficient for 
the method we propose. 

Once these calibrations are done, a number of rich clusters of galaxies at redshifts 
Zd ^ 0.3 should be chosen. Each cluster field should be divided into a number of patches, 
where each patch is large enough to contain many galaxies, but small enough that the 
lensing properties of the cluster do not vary significantly within a patch. We will quantify 
the patch size later. Within each patch, the distorted images of the faint blue galaxies 
should be used to determine the local mean values of ^g^c, Xc and (fc as functions of /. 
By comparing the local 6s,c{I) with the 6s,i{I) of field galaxies, the magnification of the 
galaxies /i(/) can be inferred as a function of surface brightness /. Similarly, from the local 
values of Xc(-^), the shear components of the lens can be determined as functions of /. As 
we explain below, this information then allows us to derive both the redshift distribution 
of the galaxies and the mass density of the lens. 

BASICS OF LENSING 

We briefly review here the relevant properties of gravitational lensing (see Blandford 
& Narayan 1992, Schneider, Ehlers & Falco 1992 for more details.) The local properties 
of a lens are described by the convergence «:, the shear 7, and the position angle (p^ of the 
shear, or equivalently by «:, 71 = 7(cos2(^-y), and 72 = 7sin(2(^-y). The convergence k, is 
deflned by 

E c^ D 

'«(^d,^s) = ^^7— — T , S„(^d,^s) = — ^ ' , (2) 

and is the surface mass density E of the lens expressed in units of the critical density Ecr- 
The quantities Dd,s,ds are angular-diameter distances from the observer to the lens and 
to the source, and from the lens to the source, respectively. Let us write the convergence 
n{z^,Zs) and the shear components ^iiz^^z^) as 

k{za,Zs) = Koo{za) f{zA,Zs) , ^t{zA,Zs) = ^i,oo{za) f{zA,Zs) , f{zd,Zs) = —-^, (3) 
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where K,oo{zd) and 700 (-^d) refer to a source at inflnity. The advantage of this split is that 
the properties of the lens are described just by K,oo{zd) and 7oo(-2d), while the variations 
induced by changes in the source redshift are described by the factor /(zd^Zs). Fig. 1 
illustrates the behavior of f{zd^ ^;s) as a function of Zg for four values of Zd- 

Given measurements of 6s,i{I) and 6s,c{I) in a given patch at a given surface brightness 
/, the magniflcation /i(/) can be calculated by 

^^^ [^"s,f(/)J ' ^ ^ 
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Fig. 1 — The function f{zd^ Zg) in eq. (3) for 
four different values of z^^ as indicated. 



while /i(/) itself is related to the corresponding k,(I) and 7(/) by 
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The ellipticity parameters xi , X2 of a galaxy image are related to the intrinsic elliptic- 
ity of the source via a transformation which depends on k, and 7. Since the orientation of 
the intrinsic ellipticities of the sources is random, it is possible to calculate from the mea- 
sured xi, X2 values of galaxies of a given surface brightness, two quantities, <^i(/), ^2{I)^ 
which describe the local distorting effect of the lens corresponding to that / (Kochanek 
1990, Miralda-Escude 1991, Schneider & Seitz 1994). The quantities <^i(/), ^2{I) are re- 
lated to n{I) and 7(/) via 
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2[1 - k{I)]^i,2{I) 
[1 - «:(/p + 7(/)2 • 



(6) 



Equations (5) and (6) can in general be inverted to estimate the convergence n{I) 
and shear 7(/) as functions of / in each patch of a lens field. In the limit of weak lensing 
(/€,7 ^ 1), these equations simplify to 



/i(/) = 1 + 2«:(/), ,5i,2(/) = 271,2(1). 
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Therefore, in this limit, the change of the angular sizes of sources relative to an unlensed 
field directly measures the convergence, while the induced ellipticity in the images measures 
the shear. 

DESCRIPTION OF THE METHOD 

Let us bin the galaxies into a number of surface brightness bins, li. For each patch in 
the lens field, and each brightness bin J^, we obtain via eqs. (5) and (6) the convergence 
n{Ii) and the shear components, 7i(/i), 72(/i). Compare now the results in two separate 
brightness bins, li and /j, within a given patch. From eq. (3) we see that 



ili) ^,{U) 72(/«) fWMU)] ~ ''^^^' 



(8) 



The quantity Rij(zd) is a kind of redshift ladder, which for a fixed z^ is a function only of 
the mean redshifts, z^ili) and Zs{Ij)^ in the two brightness bins, li and Ij. Equation (8) 
shows that, in each patch of the lens field, we have three independent estimates of Rij. 
Moreover, Rij should be the same for all patches in the lens. Thus, from iVpatch patches, 
we obtain SiVpatch independent determinations of Rij. 

Let us assume that we know the mean redshift Zs{Ii) of the galaxies in the first (i.e. 
brightest) surface brightness bin. Then, assuming a cosmographic model (e.g. Oq = 1, A = 
0), we can directly obtain from Rij(zd) the mean redshift z^ili) of all the other brightness 
bins. In any given patch the accuracy will be very poor. However, we can average over 
all the patches in a lens and also average over different clusters (which can be at different 
redshifts). Thus, the accuracy increases as more data become available. The power of the 
method lies in the fact that every new measurement can be combined with all previous 
measurements to improve the overall accuracy of the redshift determinations. 

Finally, once the redshift ladder has been used to obtain all the z^ili), then eq. (3) 
gives the lens convergence k^o and shear components 71,00, 72,00 in each patch of each lens. 
From this surface density maps of all the lensing clusters are obtained directly (e.g. eq. 
2). 

ESTIMATE OF SIGNAL-TO-NOISE RATIO 

In order to estimate the signal-to-noise ratio that we can expect with the proposed 
method, we make the following simplifying assumptions. (1) We take the number density 
of faint blue galaxies to be ~ 3 X 10^ deg~ = 83 arcmin" , which corresponds to a limiting 
magnitude of i? ~ 27 (Koo 1986, Guhathakurta, Tyson & Majewski 1990, Gardner, Cowie 
& Wainscoat 1993, Tyson 1994). We assume that the area (Ai?)^ of the observed patch is 
chosen such that the number of galaxies in the patch, ~ 83 Ai?^, is much larger than unity. 
We also assume that we have 10 surface brightness bins so that the number of galaxies per 
patch per bin is iV ~ 8.3 Ai?^. (2) We take the intrinsic scatter in the scale lengths 6^ of 
galaxies with a given / to be ±20% (Fig. 5 in Tyson 1994). The local n{I) corresponding 
to a given bin of / in a selected patch in the lens field can therefore be determined to an 
accuracy of ~ 0.2/v iV. (3) We take the rms intrinsic ellipticity of the faint blue galaxies to 
be Y (x^) — 0.3 (Tyson & Seitzer 1988, Miralda-Escude 1991). Since in the weak-lensing 
regime (x) ~ (5 ~ 27, we assume an uncertainty in the local determination of shear of 
~ 0.15/viV. (4) We model a "generic" cluster as an isothermal sphere with an Einstein 
radius of ~ 30" for a source at infinity. Thus, Koo = 7oo = 0.25/??, where ^ is the angular 
separation in arc minutes from the cluster center. (5) We assume that the faint blue 
galaxies are distributed in the redshift range z^ — 1 — 'i- If we take a cluster redshift of 
z^ ~ 0.5, the function f{z^, z^) varies by about a factor of two across this range of z^ (see 
Fig. 1). With 10 bins, the fractional change in n and 7 from one bin to the next is then 
(A/€//€oo) ^ (A7/700) ^ 0.05. 

According to assumptions (2) and (3), the noise in the determinations of n{Ii) and 
7(/i) in each brightness bin of a patch are A/€ ~ (0.2/\/8.3)A7? ~ 0.069/A7? and A7 ~ 
0.052/A7?. By assumptions (4) and (5), the actual change in n and 7 from one bin to 
the next is ~ 0.05/€oo — 0.013/??. Thus the effective signal-to-noise ratio with which the 



differences in n and 7 between two neighboring bins in a single patch can be determined is 



S_\ ^ f 0.013/^ y f 0.013/^ y ^ A^ 

To determine the redshift distribution Zs{T) of the galaxies, we average over all the 
patches in a single lens as well as over different lenses. Let us therefore integrate {S/NY 
over an area extending from the critical radius at i? ~ 0.5 out to i? ~ 5, and let us assume 
that we have data on Nd clusters. The net signal-to-noise ratio is then 



With data on say 10 clusters, we see that we can achieve a respectable signal-to-noise ratio 
of (S/N) ~ 3.5. What this means is that we will be able to determine the change in K,(Ii) 
and 7(/i) from one brightness bin to the next, and therefore the change in mean redshift 
Zs{Ii)^ with an error of only ~ 30% of the difference. Of course, if we look at the accuracy 
across the entire 10 bins, i.e. from the brightest to the faintest surface brightness galaxies, 
the fractional error will be very much smaller. 

Consider next the accuracy that we can expect in the reconstruction of the lens surface 
mass density E. Now we need to consider each patch of each lens separately since we do 
not have the option of averaging over different patches. It is therefore critical to choose the 
patch size appropriately. Assuming we have determined Zs{I) of all the brightness bins, in 
this stage of the calculations we can use all the galaxies in a given patch to determine the 
local lens parameters, Kqo, 7i,oo, 72, oo- By assumptions (2) and (3), the accuracies with 
which the local k^o and 7^0 of a patch can be determined are ~ 0.022/A7? and ~ O.OIG/Ai?, 
respectively. Thus, the effective signal-to-noise is 

5\ KqqA^ ^ Ai? fS\ 7ooA^ ^ Ai? , , 

— ^ ^ 11 — ^ ^ lo (111 

nJ^~ 0.022 ~ I? ' \nJ^~ 0.016 ~ ^ ' ^ ^ 

In the simplest scheme, we would just take the derived k^o in each patch and calculate 
the local surface density E through eq. (2). In this case, the signal-to-noise ratio is just 
llAi?/??. However, we note that k^o and 7^0 are not independent, but are related through 
a common scalar potential (e.g. Kaiser & Squires 1993). Therefore, a more sophisticated 
algorithm would combine the estimates of k^o and 7^0 across the whole lens field and 
would thereby obtain a global reconstruction of the surface density distribution. We do 
not describe the details of such an algorithm here, but note that the effective signal-to- 
noise ratio in any such scheme is likely to combine (in quadrature) the individual S/N in 
K, and 7. Thus, the (S/N) we can expect in the reconstructed E is 

S\ A{) 



This shows that if our goal is to achieve a certain target signal-to-noise ratio (S'/iV)target 
for the reconstruction of E, then we must select the angular scale of the patches to be 



A^ ^ (5/iV)target 

7? ~ 19 



(13) 



The angular resolution that we can achieve is maximum near the center of the lens and 
decreases farther out. 

Note that all the estimates in this section have been derived under the assumption 
of weak lensing. In the strong-lensing regime, especially close to critical curves, the real 
accuracy that can be achieved is likely to be better than our estimate because of the 
strong dependence of the location of the critical curve on Zg. Also, we have assumed that 
the lens is a smooth singular isothermal sphere. Actual clusters probably have significant 
substructure and this will lead to a stronger signal and therefore an improved signal-to- 
noise ratio. 

DISCUSSION 

We propose to call the redshift ladder technique described in this paper the lens parallax 
method. Although the technique does not employ true trigonometric parallaxes, it is 
nevertheless a purely geometrical method. In essence, the lens produces a differential 
bending of light rays described by the three local parameters, Kqo, 7i,oo and 72,00- The ray 
deflections distort the images of background sources in such a way that the magnitude of 
the distortions depends through simple geometry on the distance of the source behind the 
lens (eq. 3). We think the analogy with trigonometric parallax is rather close. 

There are two key ideas in this paper. First, we show how the lens parallax method 
can be used to determine the redshifts of faint galaxies as a function of surface brightness. 
Second, we suggest that by using fleld galaxies as calibrators we can measure not only the 
local shear due to a lens, as has been considered in the past, but also the local magniflcation. 
These two ideas work best when combined, but they are essentially independent of each 
other. 

The lens parallax method, for instance, can in principle be applied using only infor- 
mation on shear. This is because the quantity Rij in eq. (8) can be estimated from either 
the convergence or the shear. The signal-to-noise will of course improve if both k, and 7 
are used, but in principle 7 alone will suffice. 

We must emphasize that the lens parallax method does not make any assumptions 
about the properties of distant galaxies or their evolution. All we require is that there 
should be a reasonably tight relation between surface brightness and redshift, for which 
there is already strong evidence (Tyson 1994), and that the galaxy populations in different 
directions (i.e. towards different cluster flelds or blank flelds) are essentially the same. 
Thus the method is quite model-independent. 

The idea of using fleld galaxies as calibrators to determine the local magniflcation /i of 
the lens has been proposed in a different form by Broadhurst et al. (1994). They suggest 
calibrating the galaxy counts in blank flelds and determining how the counts vary in a 
lensed fleld. We believe our method is simpler and less model-dependent. We also suspect 
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that the signal-to-noise ratio achievable with our method is superior but this needs to be 
tested with simulations. 

As emphasized by Broadhurst et al. (1994), a measurement of the convergence n 
directly gives the local surface density E of the lens, whereas in the methods developed so 
far which are based solely on shear (Kaiser & Squires 1993, Seitz & Schneider 1994) one 
needs to convert a map of 7 into a map of E via an integral transformation. The transform 
itself is fairly complicated, but in addition the final answer is uncertain up to an overall 
additive uniform mass density (Schneider & Seitz 1994). Our technique is more direct and 
does not suffer from the mass ambiguity. 

The methods described here can be extended in several ways. First, instead of using 
only surface brightness to characterize the galaxy population, we could also include other 
properties that are invariant to lensing, e.g. colors (assuming there is no reddening due 
to the lenses). Suppose for instance that the dispersion of redshifts per bin becomes 
smaller when the galaxies are characterized by both surface brightness and color rather 
than brightness alone. It may then be possible to obtain the mean redshift of the galaxies 
as a function of both brightness and color wihtout much loss in the signal-to-noise ratio. 
We will thus achieve a more detailed description of the evolution of the faint blue galaxies. 
Another possibility is that by comparing the results from lenses at different redshifts z^, we 
can test for the internal consistency of the independent redshift ladders obtained from each 
cluster. In principle, this will allow us to constrain the model of the universe. Gravitational 
lenses are particularly effective at distinguishing models with a large cosmological constant 
A (Turner 1990, Kochanek 1993). The lens parallax method may therefore provide a purely 
geometrical technique for constraining the value of A. Finally, we can test whether or not 
the entire lens effect is due to the lensing cluster as assumed in our analysis. Under this 
assumption, the shear angle (fi'y{I) in a given patch of the lens should be independent of /. 
If the data do not show this, then we could infer that there are significant additional shear 
contributions at redshifts z > z^ which introduce differential effects as a function of source 
redshift. This could potentially be a powerful diagnostic of structure at high redshift. 

The role of systematic effects in the lens parallax method is unclear at the moment. 
The method requires that we be able to measure the surface brightnesses, scale lengths 
and ellipticities of very faint galaxies with reasonable accuracy. The effect of atmospheric 
seeing may be particularly serious, and it is possible that the method may require space- 
based observations in order to avoid errors due to variable seeing. One point to note, 
however, is that a wrong estimate of the seeing leads to errors of opposite signs in k, and 
7. It is thus possible that the method may be fairly insensitive to moderate errors in the 
seeing estimate. 

This work was supported in part by NSF grant AST 9109525. 
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